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Abstract. New feasibity of coherent quantum control of four- wave mixing processes 
in a resonant Doppler-broadened medium are studied. We propose a technique which 
enables one to enhance the quantum efficiency of nonlinear optical conversion. At the 
same time, it allows one to decrease the required intensities of the fundamental beams 
compared to those necessary in the approach based on coherent population trapping. 
The major outcomes of the analysis are illustrated with numerical simulation addressed 
within a practical medium. 

PACS: 42.50.Ct, 42.50.Gy, 42.50.Hz, 42.65.Ky, 42.65.Dr 



1 Introduction 



The concept of quantum coherence and interference plays an important role in 
resonant nonlinear optics and spectroscopy Q . Coherent coupling of laser radia- 
tion with multilevel quantum schemes has been used to manipulate energy level 
populations, nonlinear-optical response, refraction and absorption of a resonance 
medium Q. Much attention has been shown to four- wave mixing (FWM) pro- 
cesses with maximal coherence on a Raman transition [|| based on the effect of 
coherence population trapping (CPT) H|, which allows a dramatic increase of 
conversion efficiency. In far-from-degenerate schemes, substantial Doppler broad- 
ening of the resonant quantum transitions gives rise to inhomogcncous coupling 
of the driving fields with the atoms in different velocity intervals. Therefore, 
under relatively low intensities, only a small fraction of the atoms can be con- 
currently coupled with both the fields driving the Raman transition and, as a 
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consequence, the CPT conditions can not be fulfilled for the entire thermal ve- 
locity interval. This fundamental limitation in resonant nonlinear optics of gases 
substantially reduces the efficiency of the corresponding processes. The common 
methods of Doppler-free (DF) nonlinear spectroscopy entail ladder energy level 
schemes and equal frequency counter-propagating beams. This usually leads to 
large detunings from the intermediate resonance and consequently a decrease 
of the nonlinear susceptibility. Besides that, the method can not be applied to 
FWM processes in ladder schemes, due to the phase-matching requirements, and 
in Raman- type schemes due to the difference in frequencies of the coupled fields. 

The technique, allowing one to overcome this obstacle by compensating for 
Doppler frequency shifts with light-induced ac-Stark shifts and thus enabling 
one to couple a wide velocity interval concurrently, has been proposed in 
This was developed in further publications Q and most recently in publications 
[0 addressing FWM under conditions where CPT is not possible. Related effects 
in the absorption index were also investigated recently Q . 

We further mean CPT as the process where two driving fields form nearly 
equal velocity-averaged populations of the levels at Raman transitions, while the 
velocity-averaged population of the intermediate state is much less. Along with 
this absorption and amplification indices for the couple fields reach minimum. 
On the other hand coherence in the Raman transition, which is origin of FWM, 
approaches maximum. That allows one to accomplish a maximum of the con- 
version efficiency. The conditions to achieve CPT may differ from those required 
for the DF Raman-type resonance. This present paper is aimed at investigating 
what are the most favorable conditions to accomplish effective FWM frequency 
conversion of radiation at a low fundamental radiation level. This also addresses 
the problem "nonlinear optics at the level a few photons per atom" B. 



2 Theoretical 



2.1 Atomic coherence and 








Fig. 1. Energy level configuration 

ence) can be found in the same form 
driving radiations. Then the equations 



level populations 

Consider a Raman-type scheme of 
FWM, u>s = wi — u>2 + depicted in 
Fig. [I], where E\ and E2 with frequen- 
cies coi and L02 are traveling waves of 
the driving fields. The field E3 with 
the frequency L03 and generated radi- 
ation Es with the frequency cos are 
assumed weak (i.e., not changing the 
level populations). All waves are co- 
propagated. Initially, only the lower 
level is populated. The steady-state 
solution for the off-diagonal elements 
of the density matrix (atomic cohcr- 
of traveling waves as the the resonant 
or the density matrix amplitudes to the 
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lowest order in £"3,5 can be written as 

foi^oi = -i [G i(ri - r ) - G 2 ir 02 ] , P 2 ir 2 i = -i [G 2i (ri - r 2 ) - G 01 r 2 o] , 

Po2?"02 = -i [Gair 12 - Gi 2 r i] ; (1) 
An = 2Re [iGi 2 r 2 i + iGi r i], r 2 r 2 = -2Re [iGi 2 r 2 i] + 72^1, 

r + r 1 +r 2 = l; (2) 
^23^3 = -i [-G 23 r 2 + G 21 r 13 - G 03 r 2 o] , 
Po3?"03 = —i [— G 03 ro + G iri 3 - G 23 r 2] , 
Pi3ri3 = -i [Gi 2 r 23 - G 2 3?-i2 + G w r 3 - G 03 ri ] (3) 

Here and are the amplitudes of the off-diagonal and diagonal elements of 
the density matrix, Py = Py + i(Qi — kiv), fl\ is the resonance detuning for the 
corresponding resonant field (e.g., Q\ = u>i — u>oi), Pjj arc homogeneous half- 
widths of the transitions (coherence relaxation rates), v is the projection of the 
atom velocity on the direction of the wave vectors k;, 7i >2 are relaxation rates 
from level 1 to and from 1 to 2, accordingly, I\ , r 2 are relaxation rates of the 
populations of the levels 1 and 2 accordingly, G^- = —Ejdij /2h are the Rabi 
frequencies, and dij are the electric dipole momenta of the transitions. Only the 
resonant couplings are accounted for. 

The solution for the coherence r$ 2 and for those at the allowed transitions, 
which determine absorption and refraction at the frequencies u)\, lo 2 , LJ3, are 
found in the form 

_ G iG 12 f Ar w Ar 21 \ _ G m G 12 

■ Got f A |G 12 | 2 (Ar w Ar 2l \\ ,G i „ 
r oi = -i-pr- < Ar 1Q - - — — \ = -1— — 



Poi I P02 V Poi P21 J J Poi 

.G 12 J. |G 01 | 2 (Ar w Ar 21 \\ ,G 12 



P£i I P02 V Poi P 2 \ ) J P 2 *i 

.G 2 3 f / I G*i2 1 2 \ |GoiGi 2 | 2 P02 [C12 1 2 „ 
r 2 3 = »-=- m l- - — -h — „ „. - — - Ri 



P23 { v P23P3/ P03A3 PoiPii p^Pis 

where Arij are intensity-dependent population differences, and Pq 2 and Pi 3 are 
two-photon denominators "dressed" by the driving fields: 

p _ p , l G 0l| 2 , l G 12| 2 p _ p l G 0l| 2 , l G 12| 2 a s 

P)2 - P)2 H 5; 1 5 , Pl3 — Pl3 H 5 1 5 ■ (4j 

^21 P^l "03 -^23 

The coherence ro3 consists of two terms. One, fo3, determines absorption 
and refraction at the frequency u>s, and the other, fo3, determines FWM at 

UJl — 0J 2 + L0 3 — Ll>S- 

.G03 f ( A |Gqi| 2 \ |GqiGi 2 | 2 Rq 2 t |Gqi' 2 



r 3 = t— <r 1 - — -g — + P x ^, (5) 

Pl3 I V P03Pl3/ Pi3Pl3 P 0lP21 P *iPl3 



GqiGi 2 G 2 3 J / |Gqi| \ PqiP 2 *i / r 2 P : 



r 03 = ~» -p p , p i ^?02 — -g- -11+ "p- + 75T f ' (6) 

P0lP 2 iP03 L \P03Pi3 / Pl3 V-^23 P21/J 
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where P 03 = r 03 + i(Qx — i? 2 + ^3 — k s v) and k s — k\ - fe + fe. 

Making use of the above equations, the solution for the populations can be 
found from (g) as 

r2 = Y lYi -Y 3 Y 2 > ri= Y lYi -Y 3 Y 2 ' ^, = 1-^-^, (7) 



where 



and 



Fx = 2^1-^02+71, F 2 = VF 02 + VF i+r 2 , 
Is = W02- W21-72, n = W 02 + W 31 + A, 



"01 "02 
1 - "S^T" ) 



^2*1 ^02 

W Q2 = 2|G i| 2 |G 12 | 2 Rc {1/^01^2*1^02} • 



2.2 Coherence-induced Doppler-free resonance 

The appearance of a DF resonance and therefore the coupling of molecules 
from a wide velocity interval can be understood as follows. Ac-Stark shifts 
Im {IGoil 2 /^} and Im { IG12 | 2 /-Foi } hi originating from the coherence 
7"02j depend on the radiation intensity and frequency detunings. The later ones, 
in turn, depend on Doppler shifts. This allows one, by making judicious choice 
of the intensities and detunings of the driving fields, to "draw" into the dressed 
two-photon resonance all the molecules, independent of their velocities. In the 
limiting case, when the detuning from the intermediate resonance is much greater 
than the Doppler HWHM of the allowed optical transitions, the modified two- 
photon resonance described by the denominator P02 can be presented in the 
lowest order in fe, 2 t>/^i. 2 as 

P 02 = f Q2 + 1&02 fc l - + ^5^) *&} V, (8) 

. f> p l^oil 2 n |Gi 2 | 2 =l |G i| 2 |Gi 2 | 2 
where r 02 = r 2 H t^— A2 H j^-"oi, ^02 = S2i — !2 2 -\ — 

J/ 2 *'i *'2 ''l 

are the power broadened HWHM and ac-Stark shifted two-photon resonance 
detuning. From (||) it follows that the requirements for the induced Doppler-free 
resonance to be achieved arc i?o 2 = and fe = fe, where 

k\v = I 1 + — -j-^ — ) k\v and few = ( 1 + ^ ^} I feu 



are intensity-dependent Doppler shifts. Eventually the equations take form 

12 



^1 / V « 



ft(i-^-«.fi-^)=o. (10) 



/? 2 y ' v a 
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From (p|)-(^0|) follows a cubic equation as fl 2 : 



Ql r -(2-K)Ql r n 1 -Q 2r \G m \ 2 +KQ 1 \G 01 \ 2 = 7 K = k 2 /k u (11) 



which determines the detuning of field E 2 , corresponding to the induced reso- 
nance, under given values of the Rabi frequency and detuning of field E\. Here, 
the Rabi frequency of E 2 must be 



One of the roots of the equation ( |TT| ) corresponds to the dressed DF two-photon 
resonance, 



where A = 9/?i|G i| 2 (l - 2K) + 12? (8 - 12K + 6K 2 - K 3 ) + 3[3J? 2 |G i| 4 x 
x(-l - 8K + UK 2 ) + 3f2f\G 01 \ 2 (~8K+12K 2 - 6K 3 + K A ) - 3 1 G i | 6 ] 1/2 , 
B= (\G 01 \ 2 + 4f2 2 -A[2 2 K + Q 2 K 2 )/'dA 1 / 3 , K = k 2 /k 1 . 



The conditions for the CPT and DF resonances may fit each other or differ, de- 
pending on the specific case. This and corresponding outcomes regarding FWM 
will be investigated below numerically for the most optimum situations, where 
the analytical solution can not be obtained. 

3 Numerical analysis 

The graphs presented below are computed numerically based on averages over 
the velocity equations (||) , (||) and (^) . For the numerical analysis we have used 
the parameters of the sodium dimer transition with the following wavelengths 
|0): X 1 S+(v" = 3, J" = 43) - B 1 7T„(6, 43)(A i = 488 nm) - X 1 £+(13,43) 
(A 2 i = 525 nm) - A 1 S+{2A, 44)(A 23 = 655 nm) - X 1 r+(3,43)(Ai 3 = 598 nm). 
The corresponding homogeneous half-widths of the transition arc 20.69, 23.08, 
18.30 and 15.92 MHz, whereas the Doppler half-widths Aujiy) are equal to 0.92, 
0.85, 0.68 and 0.75 GHz. The numerical simulations allow us to analyze the 
velocity-averaged equation for the dctunings, where the approximation taken in 
(H) is not valid. 

3.1 Nonlinear resonances in FWM polarizations, in absorption 
indices and in the level populations 

Figure |[ plots 1 illustrate the coherence-induced compensation of Doppler shifts 
with ac-Stark shifts resulting in the resonance narrowing in the squared module 
of the velocity-averaged reduced nonlinear FWM polarizations 



W = (|X4| 2 /IX40| 2 )|ff0l| 2 |5l2| 2 [Fig! ( & )] and l^l 2 = (IX3|7|X40| 2 )|<70l| 2 k2 



[Fig.g] (b)]. The nonlinear susceptibility is reduced by its maximum value in 
the same frequency range, but for ncgligibly-wcak E\ and E 2 fields, goi = 




(12) 



fl 2r = (A 1/3 + 9B + 2f2 l - QiK)/3, 



(13) 
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x 10~ 3 




Fig. 2. Normalized squared module of the FWM polarizations P4 (a) and P3 (b) 
vs the detuning of driving field E2. The plots correspond to: 1- DF resonance, 2 
- CPT, 3 - intermediate case, subplot - same resonance, but under weak fields. 

Gqi/Aujid, gi2 = G12/ 'Alu2d- A substantial narrowing is seen from comparison 
with the subplot, presenting the spectral dependence of the velocity-averaged 
squared module of the same susceptibility, but in the negligibly weak fields, nor- 
malized to unity. The HWHM of the resonance in the subplot is approximately 
70 MHz, which corresponds to the Dopplcr width of the nonperturbatcd Raman 
transition, whereas the HWHM of the resonance in plots 1 is 15 MHz. This 
indicates some power and residual Doppler broadening. The Rabi frequencies 
of the fields E\ and E 2 are equal to 157 MHz and 85 MHz. The detuning of 
Ei is equal to 413.8 MHz (0.45 • Aujid), and detuning of E 3 is equal to -695.5 
MHz (—1.02 • Aujzd)- Plots 2 display the same dependences, but related to the 
CPT conditions (the half- width of this resonance in I-P4J 2 is 200 MHz). In this 
case, the required Rabi frequencies of E\ and E2 are equal to 351 MHz and 332 
MHz, respectively, while the detunings are identical to those in plots 1. Plots 3 
display the same dependence, but under the intensities, that provide narrowing 
in the intermediate range (HWHM of resonance is 70 MHz). In this case, the 
Rabi frequencies of the fields E\ and E2 are equal to 222 MHz and 235 MHz, 
respectively, whereas the detunings arc the same as in the previous plots. The 
important outcome is that in the Doppler-free regime, larger FWM polariza- 
tions can be accomplished under lower driving fields. Also, we want to stress 
the substantial difference in magnitudes of the nonlinear susceptibilities \4 an d 
X3 , which follows from the interference of contributions of molecules in different 
velocity intervals to the macroscopic nonlinear polarizations. Figure |^(a) shows 
dependence of the velocity integrated populations on the detuning of the field 
E2, while the Rabi frequencies of E\ and E2 are equal to 351 MHz and 332 MHz, 
respectively so that the CPT conditions can be fulfilled under appropriate de- 
tuning. The subplot displays the distribution of the populations over velocities, 
whereas the detuning corresponds to the CPT regime for the velocity-integrated 
populations. Here u is the thermal velocity with Maxwell envelope removed. 
The subplot indicates CPT for molecules in v — 0; however, even inversion of 
the populations occurs in a relatively wide velocity intervals. 

Figure |^(b) is computed for the Rabi frequencies of the fields E\ and E2 
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Fig. 3. Velocity integrated population of the energy levels vs detuning of driving 
field i?2- a- near CPT and CPT regime, subplot - velocity distribution of the co- 
herently trapped populations; b- near DF regime, subplot - velocity distribution 
of populations under DF regime. 

equal to 157 MHz and 85 MHz, respectively. The inset displays the distribution 
of the populations over velocities at the detuning corresponding to DF. The 
distribution is even more complicated compared to that in the previous case. 

Figure [| displays velocity-averaged absorption indices vs detuning of the field 
Ei . The intensities of the fields are the same as in the corresponding graphs given 




Fig. 4. Velocity averaged absorption indices ai,2,s near DF (a) and CPT (b) 
regimes vs detuning of the driving field E^. 

above. The graph |] (a) shows a substantially decreased absorption index for the 
first field due to saturation effects, which is almost independent of the detuning 
of the second field because of its relative weakness. The second field experiences 
Stokes gain. The plot for the generated beam (reduced by 10) displays a nonlinear 
resonance in the absorption index. With growth of the second field intensity, so 
that CPT becomes possible, both absorption of the first and gain of the second 
driving field experience substantial decrease in the CPT regime. Absorption of 
the generated field dramatically changes as well (Figure H (a)). 

A linear approximation in ki 2^/^1,2 (see equation (|8j)) is valid for large de- 
tunings from the intermediate resonance. While for tuning close to the resonance, 
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the concurrently coupled velocity interval decreases, the nonlinear response of 
the molecules grows. At the same time, absorption of the coupled radiations 
increases as well. 

The output FWM generated radiation is determined by the interplay of the 
above considered processes. Moreover, resonant FWM may not be viewed as 
a sequence of the independent elementary acts of absorption, gain and FWM. 
Interference of these elementary processes plays a crucial role, even for a quali- 
tative prediction of the spatial dynamics of the generated field pi] . Therefore, 
appropriate optimization is required in order to achieve the maximum FWM- 
generation output. In the next section, relevant results are presented based on 
the above analyzed dependencies. 



3.2 Coherent quantum control of FWM in a double- A Doppler 
broadened medium 

The solution of the Maxwell equations for the coupled traveling waves can be 
found in the form 

E ] (z,t) = Re{E j (z)exp[i(uj j t - kjz)]}, (14) 

where kj are complex wave number at the corresponding frequencies, kj = k'- — 
iaj/2. The set of the equations for the amplitudes, relevant to the case under 
consideration, is 

-dz-~~T Eu ~dz-~-Y E2 ' (15) 

^ = - y £ 3 + a 3 E\E 2 E s ex V {-iAkz}, (16) 

^ = E s + a s E x E* 2 E 3 exp{iAkz}, (17) 

where a 3 = i2irk 3 x 3 and &S = i^ksxs are cross coupling parameters, \3 an d 
Xs are intensity-dependent nonlinear susceptibilities for the FWM processes 
us «-> toi — uj 2 + UJ3 and Ak = k$ — fci + k 2 — k% . The boundary conditions are 
Ei(z = 0) = E i (I = 1, 2, 3) and Es(z = 0) = 0. We assume that intensities of 
the weak waves with frequencies uj 3 and ojs are weak enough so that the change 
of the strong fields E\$ due to FWM conversion can be neglected. The quantum 
efficiency of conversion QEC of the radiation E 3 in Eg is defined by the equation 

Vq = (uj 3 /u; s )\E s (z)/E 3 (0)\ 2 exp(~a s z). (18) 

Assuming that the condition of phase matching Ak' = can be ensured (e.g., 
with a buffer gas), the solution for QEC of the above given equations can be 
expressed as 



uj 3 I as i2i 1 2 



r li = ~\~\ \EqiEq2\ exp (Rego--)z |1 - exp(-g ^)| , 



i2 



5o 



2 , A I T7i T7i |2 



[(a>i + Q' 2 + Q'3 - as/2)] + Aasa3\E iE { 



1)2 



1/2 



a 3 + as- 
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These formulas account for nonlinear resonances both in FWM nonlinear polar- 
izations and in absorption indices. In order to compute the absolute magnitude 
of QEC (under the assumption of Ak' = 0) we have used the data for the Frank- 
Condon factors of 0.068, 0.142, 0.02 and 0.036 for the transition with wavelengths 
Aoi, A 2 i, A23 and A03 accordingly. 

Figure |^ presents a numerical simulation of the evolution of QEC along the 
vapor cell. The distance is scaled to the resonance optical density amz with 
the driving fields being turned off. The plot 1 corresponds to the case where 
the conditions for Doppler compensation at the transition 0-2 are fulfilled on 
the entrance of the cell. The Doppler-free resonance would be narrower if one 
of the coupled fields were weak. However, that would give rise to larger absorp- 
tion and to decreased nonlinear FWM polarization. The optimization accounting 
for the change of the driving fields along the medium also shows that it does 
not substantially change the optimum input intensity values for the plots pre- 
sented in Fig. ||. Plot 2 displays the same dependence, but for the optimum CPT 
conditions. Alternatively, they are not optimum for the elimination of Doppler 
broadening. 

The maxima of these two curves are 
comparable, but in the Doppler- free 
regime, quite lower intensities of the 
fundamental radiations are required. 
Plot 3 displays a similar dependence 
in the intermediate case. All inten- 
sities and detunings are identical to 
those used for computing the previ- 
ous figures. 

For the transitions under consid- 
eration, the Rabi frequency on the 
level of 100 MHz corresponds to pow- 
ers of about 100 mW focused on 
a spot on the order of 10~ 4 cm 2 , 
which can be realized with common 
cw lasers with the confocal parame- 



5i 

r| ,% 








3 /\ 
5 / 






- 1 

2 




2 / 






3 




5- // \ 








1 - 








5-1 








0'- 









10 



15 



Fig. 5. QEC vs optical density of the 
medium. 1 - DF, 2 - CPT, 3 - interme- 
diate regimes. 

ter of focusing of about 2.5 cm. At a temperature of about 700 K, the optimum 
optical density of the dimers is achievable for vapor lengths of about 2 cm. This 
is in accord with the parameters of typical cw FWM experiments (see, e.g. fl(i[| ). 

The above indicated intensity corresponds to about 10 11 photons per cubic 
centimeter. Accounting for the molecule number density at the indicated temper- 
ature, which is on the order on 10 13 cm -3 , we have about 10 -2 photon/molecule. 



4 Conclusion 

Coherent control of populations [ p"2"| and of four-wave mixing fill with pulses 
shorter than the dephasing time T2 has proven to be a powerful tool for manip- 
ulating nonlinear-optical and chemical properties of free atoms and molecules. 
In these cases, maximum coherence can be achieved as a result of Rabi oscilla- 
tions of the two-photon transition, and the required driving intensity is much 
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higher than that proposed in this paper. Consequently, Doppler broadening of 
the coupled transitions does not play an important role. 

On the contrary, this paper considers coherent quantum control of resonant 
four-wave mixing processes, where Doppler-broadcning of a double-vl medium 
is the factor of crucial importance. An approach enabling one to enhance the 
efficiency of nonlinear optical conversion and, at the same time, to decrease the 
required intensity for the fundamental pump beams is proposed. This approach 
is based on the elimination of Doppler broadening of the resonant two-photon 
transition. The advantages of the proposed method as compared to those based 
on coherent population trapping, where inhomogeneous broadening may play 
important role too, arc illustrated with numerical simulations for realistic exper- 
imental schemes. The results obtained contribute to the field of nonlinear optics 
at the level of several photons per atom, which is currently attracting growing 
interest. 
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